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Abstract 

We present a complete 1-loop study of the process 77 — > Xi'Xj and the predicted cross 
section in a 77 Linear Collider. A suitable numerical code PLATONlc, valid for any set 
of real MSSM parameters, is released. This study and code are complementary to those 
suitable for Dark Matter detection through the inverse process x^x^ —>■ 77 describing 
neutralino-neutralino annihilation at rest, which were presented previously. If SUSY is 
realized in Nature, both codes should be very helpful in future astrophysical and collider 
studies of the neutralino sector. 

PACS numbers: 12.15.-y, 14.80.Ly 
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1 Introduction 



A tacit candidate particle for the dominant contribution to the cold Dark Matter (DM) 
that apparently constitutes almost a third of the energy density of our Universe is, of 
course, the lightest neutralino(s) predicted in an R-parity conserving minimal supersjmi- 
metric model (MSSM) [H El Ej- Provided the MSSM couplings have the appropriate 
values, the annihilation processes XiXi ~^ 77; 7^ could produce observable rates of very 
energetic photons coming from the center of our Galaxy jH 13 E] ■ 

Such observable photon rates could be realized due to a large Wino or Higgsino com- 
ponent for the lightest neutralino Xi? and possibly also from resonance effects induced 
by the or masses jH El El El ■ In any case, the observability of such halo galactic 
photons depends so strongly on the MSSM parameters, that even the non-observation of 
any signal could produce useful constraints. 

Complementary, and in principle much more detail information on the neutralino 
properties could be obtained by studying e~e~^ Xi'Xj and 77 'XaX'j at a Linear Col- 
lider (LC) [HI, for any neutralino pair = 1,...,4) accessible by the available energy. 
Since the first of these processes is reahzed already at the tree level, studying its signa- 
tures should eventually supply most of the experimentally accessible information on the 
neutralinos (TUl El . 

Nevertheless, the study of the 1-loop process 77 XiX'j in a 77 Linear Collider (LCyy) 
\1'2\ . is also useful, since it directly tests in an Earthy experiment, exactly the same process 
as in the Dark Matter (DM) searches. Compared to the DM studies, the extra advantage 
of Collider measurements though, is that they can be done for any neutralino pair - even 
unstable ones, and be performed over a considerable range of energies. Furthermore, 
depending on the polarizations of the e^-beams and the laser photon beams, six different 
" cross-section-like" observables are in principle available, even if we sum over all possible 
helicities of the final neutralinos; see e.g. ^31- This is much richer than in the case of DM 
studies where only the unpolarized total ~^ 77 cross section for the cosmologically 
stable Xi neutralinos, very close to threshold, is relevant. 

The purpose of this work is to present such a study based on the complete set of the 
contributing 1-loop diagrams. In Section 2 we discuss the general properties of the ampli- 
tudes of the processes 77 — >• XiX^ and its reverse, as well as those of the LC^^-observables, 
when the neutralino helicities are summed over. The results for these observables in an 
extensive set of MSSM benchmark models are discussed in Section 3, where the released 
PLATONlc code is also presented. Finally Section 4 contains the conclusions. 



(Xi')X^) (with i,j = 1, ...,4 being the neutralino counting indices), and (/ii,/i2) the hehc- 




We denote by K 
center-of-mass sc 




(9), the invariant amplitudes for the process Xi'Xj ^ 77 at a 
; angle 6. Here (A,-, A,) denote the helicities of the neutralinos 
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ities of the photons^. The amphtudes for the reverse process 77 — > Xi'Xj ^-^e written as 
-^Mi,/z2;Ai,A2(^)- 

Following the Jacob- Wick conventions the Xi'xrfermion antisymmetry implies^ 



while the 77-boson symmetry requires 

^,..;A„A.W = (-1)'^-'^^..,A..A.(^-^^) • (2) 

If the MSSM breaking parameters and the Higgs parameter fi are real, then time 
reversal and CP invariance hold for the neutralino processes, implying also 

^Ml,M2;Ai,A2W = ^Ai,A2;mi,A.2W " (3) 
n.-.2;-A..-A.(^) = (-l)^-^-(^-^^) V^VjFl^,,X,,,M , (4) 

where r]j = ±1 is the CP-eigenvalue of the neutralino'^ x'j- 
Combining (HHHH), we get 

F\i\2;/ii/i2^^^ ^ ^^^^^ '^^-^A2Ai;^2Mi'^^) ^ '^*^J-^-A2 -Ai;-M2 -Mi^^) ' 

^.2;A.A2 W = (-l)^^-^^^'^-'^^:..A2A. W = ^.^.n2,-..-A2,-A. W , (5) 

where the first part comes from 12)) alone, while for the last part the CP-invariance 

relation (jH) is also used. On the basis of (0) we could select F++ , 

F+_4._, ^. as a possible set of independent helicity amplitudes at each specific angle. 

Restricting to real MSSM parameters from here on, it is thus sufficient to only consider 
the amplitudes of the process Xi'x^ ~^ 77? for both DM and LC^^ studies, using the 
one-loop diagrams generically described in Fig^ in the 't-Hooft-Feynman gauge. The 
types of particles running clockwise inside each loop in Fig|T] are: box"^ (a): (fSSS), 
ifWWW), ifSSW), ifSWW), ifSWS), ifWWS); box (b): (Sfff), {Wfff); box 
(c): (SffS), (WffW), (SffW), (WffS); initial triangle^ (d): (SfS), (WfW); final 
triangle^ (e): (WWW), (SSS), {SWW), (SSW); final triangle (f): (///); and bubbles 

iThus A,;j = ±1/2 and ^1,2 ^ ±1, so that (-1)a'i-M2 ^ i 

^ These relations do not agree with the antisymmetry condition used in |15| . 

■^We follow the same notation as in e.g. |llj . 

*ln naming the particle-strings in the box-loop, we always start (moving clockwise) from the line 
ending to the x'j vertex. 

^In naming these particle-strings in Fig^, we always start from the line ending to the x^-vertex. 
^The particle-string names in Fig^,f, always start from the line leaving the vertex where the s-channel 
exchanged neutral particle Z, ^ ^ or ends. 
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(g): (SS), (WW); (h): (WS). By S we denote the scalar exchanges (Higgs, Goldstone 
and sfermions); and by / the fermionic ones (leptons, quarks and inos). Bubbles (g) and 
final triangles (e) and (f), are connected to the initial state through an intermediate 
Z, or neutral Higgs or Goldstone boson h^, H^, A^, G^. 

For book-keeping the Majorana nature of the neutralinos, we always describe by 
a positive energy Dirac wavefunction, and Xj by a negative energy one. We then note 
that the triangle and bubble diagrams in Fig^ separately satisfy the neutralino and 
photon symmetry conditions (012)). For the separate box-diagrams in FigQ though, the 
neutralino exchange symmetry (P) is only satisfied after adding to each box of Fig^ the 
corresponding photon exchanged one. The calculation of these later boxes may be avoided 
though, by imposing (j2)) to the contribution of each of the boxes in FigQ Requiring 
then the validity of for each such box contribution, provides a stringent test of the 
computation. We have checked that all these constraints are exactly satisfied. 

We have derived explicit analytic expressions for the contributions of the various di- 
agrams of Fig^in terms of Passarino-Veltman functions [16J. These expressions are so 
lengthy though, that it is useless to present them explicitly. We have therefore instead 
chosen to release a numerical code, which calculates the differential cross sections defined 
below, for various photon polarizations in any MSSM model with real parameters at the 
electroweak scale; see discussion below. 

Before turning to this though, we add two remarks concerning the behavior of the 
amplitudes near threshold, and at very high energies and angles, respectively. Very close 
to threshold, the relative orbital angular momentum / of the neutralino pair XiX^ is, of 
course, vanishing. Because of Fermi statistics, if i = j, the neutralino-pair can then only 
exist in an ^So-state corresponding to a total neutralino spin ^tot = 0; while if i ^ j, 
Stot = 1 is also possible. In both cases, the CP eigenvalues of the x'iX^-pa-iic is 

Thus, for / = 0, Stot = and r/jZ/j = 1 (r/jZ/j = —1), the neutralino pair has the same 
quantum numbers as the {H^) neutral Higgs-state, leading to a resonance enhancement 
for the Aj = Aj-amplitude, provided the appropriate neutral Higgs mass is close to the 
sum of the neutralino masses^. This is induced by the diagrams in Figs^,f,g,h. 

Next, concerning the behavior at very high energy (s ~ |t| ~ |m| ^ MsusY-s>caAe^), 
the dominant leading-log contribution to the neutralino production amplitude is found to 
be generated by the diagrams in Figs^,b,c,d when the internal loops involve chargino 
and VT-exchanges. This is given by 



— ^ —sin 



^ I — (— I I ■ I I ■ I I I - I I ■ 

+ ■■■ , (6) 



a- . „ _ 1 + ri.rij , ^ 1 - ViVj 



2 2 



''Since the neutralino masses are heavier than the Z-mass in ah contemplated MSSM models, there 
can never be an analogous enhancement due to the intermediate Z boson in the i ^ j, Stot — 1 case. 
*For such kinematical situations: t ^ — s(l — cos6')/2 and u ^ — s(l + cos6')/2. 
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where the dots stand for subleading constant and 1/s-terms, while 



Li? = (/i2 - Aii)(Ai - A2)^ + 2;Ui/i2(Ai - A2) cos6' -^(^In 
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and 



Cj = 2Um{ZZZ^n +llrn{ZZZ^: + Z^^Z§:) 



(9) 



with Z^- describing the neutrahno mixing matrix in the notation of ^Tj. As expected 
from [18J, gauge invariance ehminates all isolated single-log and log-squared terms in (jHllZl 
IHI), so that the only allowed leading-log contribution consists in differences of log-squared 
terms, i.e. a single ln(s/m^) multiplied by angular dependent but energy independent 
logarithmic coefficients for the real part of the amplitudes, and simple 7rln(s/mp^/)-terms 
for the imaginary parts. At large angles, this has the tendency to enhance the imaginary 
parts of the amplitudes, as compared to the real parts. Such an asymptotic behavior has 
also been observed (occasionally very strikingly), for the helicity amplitudes of process 
like 77 — >• 77, Z7, ZZ, y4°y4°, always being generated by W-loops ITU]. 

Note (from the above expressions of and that, in the high energy limit (when 
masses are neglected and s + 1 + u = 0) the only non vanishing helicity amplitudes are 

those with Ai = — A2 and /ii = —^21 i-e. -F+_4._ and ^ and the ones related by 

symmetry relations eq.(0|S)). This is due to the dominance of box diagrams with chirality 
conserving gauge couplings of W bosons. 

In practice these asymptotic expressions have some interest only when the energy is 
sufficiently higher than the SUSY scale; for example in the TeV range if the SUSY masses 
are of the order of a few hundred GeV. At lower energies, two types of effects modify 
the asymptotic behavior. One is purely kinematic and due to the additional constant 
and energy decreasing (M^/s) terms denoted by the dots in (jHI). The other one is the 
appearance of chirality violating amplitudes due to box and triangle diagrams involving 
scalar couplings (neutralino-sfermion-fermion and neutralino-neutralino-Higgs couplings). 
These amplitudes vanish like M' / s at high energies, but at energies not too far from 
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the SUSY scale, their relative importance, compared to the one of the leading chirality 
conserving ones, is model-dependent. In particular, it depends on the neutralino contents 
(Bino, Wino and Higgsino mixture) which controls the relative magnitude of its gauge 
and scalar couplings. 

All these properties can be checked by making a numerical comparison with the com- 
plete results. We come back to the discussion of these in the next section, where we 
discuss our results for various benchmark models. 

We next turn to the 77 — > Xi'Xj collisions in a 77 collider LCyy [I^ realized through 
backscattering of laser photons in an e~e^ LC jH]. For real MSSM parameters, where 
CP-invariance holds, the general form of the neutralino production cross section is 

cicr(77 ^ X°X°) dL^^\ da^ da22 



drd COS 9 dr {dcosO dcosO 

+(&5i>:^sin2.#.+ «20:^sin2^'| ■ (10) 

a cos a cos J 

In (fTUj) . T = s/see, with See being the square of the e~e+ center-of-mass energy and 
s = = s^o^o denoting the corresponding quantity for the produced neutralino pair. 

The quantity dL^^y/dr describes the photon-photon luminosity per unit e~e'^ flux, while 
the parameter-pairs (^2, ^2)5 (^35 ^3) (0, (p') describe respectively the average helicities, 
transverse polarizations and azimuthal angles of the two backscattered photons. These 
are in turn determined by the corresponding quantities of the laser photons, and the 
e^-polarizations [121 CH] • 

Finally, the cr^-quantities in (fTUI) are defined as 



dao 
dcos9 



( i287rg ) \^++^i>'2\^ + 1-^— AiAaT + |-^+-AiA2r + |-^-+Ai 

A1A2 



|2 1 
A2I J ) 



+ + |2 + |F+_+_|2 + |F+__+n , (11) 



da22 
dcos9 



(li^) Et l^++A.A2r + Ii^-A.A2P - \F- 

V647rs J 



+-AiA2r ~ |-^-+AiA2l 



A1A2 

-(|F+_++|2 + |F+___p + |F+_+_p + |F+__+n] , (12) 
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dcosO 



da's 
dcos9 



dcosO 



~- (^4^) X][^++AiA2^-+AiA2 +^+-AiA2^--AiA2] 
A1A2 



(q^-) Xl[^++^i^2^+-AiA2 +^-+AiA2^--Ai 
^ ' ^ A1A2 

Re[F+ + F+_+_Fl 



A2J 



V327rs 

\32TrsJ 
\32ttsJ 



++- 



Re[F+ 



-++ 



(13) 



(14) 



(15) 



da' 



33 



dcos9 



da2z 
dcosO 



da'o 



23 



dcos9 



Cij 



32ns 

32TTS 



-++ 



-F, 



^++-F;__^ - F++_+F;_^_]rj,rjj , 



(3A 

\32tts J 



Im[F++++F;_^^ + F++__F;_ 



+F. 



(16) 



(17) 



(18) 



where the reduction due to the identity of the final neutrahnos whenever i — j, is taken 
into account through the coefficient 



Cij — Pij 1 



Pij — 



rrii — mj)^] r {mi + mj)^ 



(19) 



and the cm. scattering angle 9 varies in the range —1 < cos9 < +1. 
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On the basis of (^121) we obtain that dao/d cos 6, da22/dcos6, da-^^/d cos 6 and da' 33 /d cos 6 
are symmetric under the interchange 

whereas 

d(T3 

dcosO e 

d(T23 

dcosO e 

3 Results 

On the basis of the diagrams in Fig^ we have constructed a numerical code called 
PLATONlc which calculates the differential cross sections in (jlllll8|) for any cm. energy 
and any set of real MSSM parameters at the electroweak scale, using [20] • This code, 
together with an explanatory Readme file, can be downloaded from pi] . 

We have made a run for the 31 benchmark models already used in the previous DM 
study in [7j. For each model we have computed both the angular distribution and the 
energy dependence of the various cross sections presented in Section 2, for x\x\i X1X2 ^^"^ 
X2X2 production. 

The typical feature of the angular distribution is an isotropy at low energy, and a 
tendency to present some forward-backward peaking at high energies (see the 1/t and 1/u 
terms in eq. ((71 |H1) ) • For the integrated cross sections, one finds energy dependent structures 
at s ~ (mi + m2)^, with mi 2 being the masses of intermediate particles contributing to 
the diagrams of Fig^ as well as resonance effects due to diagrams involving AP or 
intermediate neutral Higgs-states. 

These features are common to several benchmark models, so we will not illustrate all 
of them. Moreover, since the expected luminosities at the future colliders will probably 
be of the the order of 10^ fb~^ per year, we ignore those model-cases considered in Tj, 
where the cross sections are much smaller than ~ 0.1 fb. We thus present results for 
the universal mSUGRA benchmark models [22] SPSlai, SPSA [22], AD{fg5)i [211; the 
non-universal mSUGRA models CDGj^ and CDGqii [23; and the GMSB models [20] 
SPS7i, SPS81 [23j. The grand scale defining parameters for all these cases are shown in 
Table 1. 

Starting from the grand-scale parameters in Table 1, their electroweak scale values and 
the and widths are calculated using the codes SuSpect [27j and HDECAY [2H1- In 
turn, these values constitute the input needed for the PLATONlc code [21] • The obtained 
results are indicated in Figs|21IHlbelow. In all cases, the sub-figures labelled (a), (b) and (c) 
describe the integrated cross sections of Eqs. (fTTllTH|) . in the range (1° < 6^ < 179°) for XiXi; 
X1X2 X2X2 production respectively. Of course, only X1X2 X2X2 production are in 
principle observable, since the lightest supersymmetric particle LSP Xi should be invisible 
in an R-conserving theory. In addition, sub-figures (d) and (e) give the unpolarized 
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Table 1: Input parameters for the used MSSM benchmark models at the grand scale. 
Dimensions in GeV. In all cases n > 0. By convention, M2 > is used. 





universal mSUGRA 


non-universal mSUGRA 


GMSB 




SPSlai 




SPSA 


CDG75 


CDGoii 




SPS7i 


SPS81 


Ml 


250 


400 


300 


-400 


424 


^mess 


80000 


120000 


M2 


250 


400 


300 


240 


200 


MsUSY 


40000 


60000 


Ms 


250 


400 


300 


80 


40 


tan/3 


15 


15 


mo 


100 


220 


400 


1400 


1400 










100 


220 


400 


1400 


1400 








Ao 


-100 








1000 


1000 








tan P 


5 


40 


50 


50 


50 









differential cross section dao/dcosO at cm energies 0.4 and 2 TeV respectively, for the 
X^x'j production channels which are energetically accessible. On the basis of Figs|21IHl we 
remark the following: 

In the universal mSUGRA cases of FigsI2E]and the non-universal mSUGRA CDG-j^ 
of Fig|31 only the ctq, cr22 and occasionally the (73 = 0-3 cross sections for 77 — ^ X2X2' 
sufficiently large to be (in principle) observable in a conceivable future Collider. This is 
because X2 ~ VF*-^^ and x? ~ -B, in these models. The largest and easier to measure cross 
section is, of course, the unpolarized ctq, defined in ffTTj) . As seen from Figs|21infc) though, 
in an LC^^ collider at 0.4 TeV, the X2X2"P^ir production can only occur in the case of 
the SPSlai and CDGj^ models. In all other cases, X2 too heavy to be produced at 
0.4 TeV, but it can be generated in a 2TeV 77 version of CLIC; compare FigsEJE^d). As 
seen there, dao/dcosd is very fiat at 0.4TeV, and develops a moderate forward-backward 
peak already at 2TeV. 

Figini concerns the OH non-universal mSUGRA model of [23] , where Xi is predom- 
inantly a Wino with some appreciable Higgsino components, while X2 is predominantly 
a Bino. Depending on the LC^^ energy and fiux, both X1X2 X2X2 productions may 
be observable. A special feature of this model is that X1X2 production is predicted to be 
generally more copious than the X2X2 o^^- Although, dao/dcos6 is very fiat at 0.4TeV for 
all channels, it develops strong forward-backward peaks by the time the energy reaches 
2TeV; compare subfigures (c) and (d). 

Concerning the GMSB models j2n] in Figs ]7l81 we remark that their general structure 
is reminiscent of the results in the mSUGRA models of Figs|21lll Finally, no results are 
presented for the AMSB models [5^ I23j . since the almost exact relations Xi ~ ^^^id 
X2 ~ -B there, enforce unobservably small values for the X1X2 X2X2 production cross 
sections. 

We next turn to the comparison of the numerical results from the exact 1-loop compu- 
tations, with those from the asymptotic expressions in ()6I71 IH)) . First, we observed that in 
all benchmark models considered above, the chirality violating amplitude and the 
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chirality conserving and |_, are the dominant ones at energies above ITeV. 

In all cases, -F_^+++ is strongly decreasing with energy in the several TeV range; while 
I and |_ become increasingly dominant (with increasing energy) approach- 
ing their asymptotic expressions from (0), up to a constant. Since (0) is generated by 
chargino-W loops, the accuracy of the asymptotic expressions at high energies is particu- 
larly striking for cases where the produced neutralino is mainly a Wino. A best example 
of these is illustrated in FiglHlfor X2X2 production in the AD{fg5)i model. There exist 
cases though, particularly in the production of Bino-like neutralinos, where the constant 
departure between exact and asymptotic expressions for the and ^. ampli- 
tudes, and the subleading amplitude, are still important at the few TeV energy 
region. 

In practice it means that an amplitude analysis of experimental data at high energy 
could immediately determine the nature of the produced neutralinos. If a precise ampli- 
tude analysis is not possible, at least a fit of the unpolarized cross section of the type 
a + 61n(s/m^) + cln^(s/m^) should produce a measurement of the logarithmic slope c, 
which could then be compared to the model predictions for the square of the quantities 
Cji or C/ in 0, thereby providing clean tests of the neutralino properties. 

4 Conclusions 

In this paper we have presented our work on the properties of the process 77 'Xa'Xj, 
that are measurable in a 77 linear collider LCy^. It complements the study of the reverse 
process in which was adjusted to the environment of the Dark Matter searches. In 
both cases, the complete set of Feynman diagrams has been calculated. 

As in jTj , we have found that these processes are very sensitive to the actual values of 
the various MSSM parameters, and details of the neutralino mixings in particular. The 
basic size of the cross sections can vary by several orders of magnitude, depending on the 
contents of the produced neutralinos, the highest values being obtained when neutralinos 
are Wino-like. Depending on the SUSY model, these properties reflect in a variable way in 
the Xi X2 X2 X2 productions. Further enhancements can appear locally owing to the 
occurrence of threshold effects due to relatively light intermediate states in box diagrams 
or of A^, resonance effects. 

In the low and medium energy range these processes depend on most of the MSSM 
parameters (because of the variety of particles appearing as intermediate states in the 
one loop diagrams of Fig.l) but we have also shown that, in the high energy limit, the 
amplitudes tend to a simple logarithmic form with a coefficient that only depends on the 
neutralino mixing matrices without any other SUSY parameter. 

The angular distribution of the differential cross sections is generally rather fiat at low 
energy, but a simple forward-backward peaking is progressively generated according with 
the asymptotic rule mentioned above. 

These processes have a very rich supersymmetric contents (in fact all SUSY particles, 
except gluinos, contribute) so that measurements at variable energy should provide very 
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stringent tests of MSSM and the neutralino structure. 

Concerning in particular the Linear CoUider measurements, we should remember that 
the LC^^ study of the 77 x?Xj-process, complements the study of the tree level LC 
observables for e~e^ X^X^j- If MSSM is really vahd, and neutralinos turn out to be 
accessible to future LC and LC^^, the combined analyses of such experiments can increase 
our knowledge on the nature of neutralinos and of SUSY in general. 

But the most important fact is that they can be combined with the dark matter 
searches, thus allowing exchange of knowledge between accelerator particle physics and 
astrophysics. To facilitate future work on LC or DM studies, the numerical codes PLA- 
TONlc and PLATONdml have been released, calculating the Linear Collider differential 
cross sections and the Dark Matter rates respectively 21]. These codes are applicable to 
any MSSM model with real parameters at the electroweak scale^. 
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Figure 1: Feynman diagrams for XiXj 77- internal lines denote fermionic ex- 
changes; while broken internal lines denote either scalar or gauge exchanges, except in 
the diagram (h), where the W and Goldstone exchanges are indicated explicitly. 
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Figure 6: Same caption as in Fig|21for the benchmark model CDGqii |25j . 
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Figure 7: Same caption as in Fig|21for the benchmark model SPS7i [221 ■ 
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Figure 8: Same caption as in Fig|21for the benchmark model SPS81 
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